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Abstract 

Recently the AMS-02 experiment reported an excess of cosmic ray antiprotons 
over the expected astrophysical background. We interpret the excess as a signal 
from annihilating or decaying dark matter and find that the observed spectrum is 
well fitted by adding contributions from the annihilation or decay of dark matter 
with mass of O(TeV) or larger. Interestingly, Wino dark matter with mass of 
around 3 TeV, whose thermal relic abundance is consistent with present dark matter 
abundance, can explain the antiproton excess. We also discuss the implications for 
the decaying gravitino dark matter with R-parity violation. 


1 Introduction 


The existence of the dark matter (DM) has been conhrmed by various cosmological ob¬ 
servations [1], yet its identity is a complete mystery. The DM provides the most robust 
evidence for physics beyond the Standard Model. 

Recently, the AMS-02 collaboration has reported their latest results of the cosmic-ray 
antiproton measurement [2], which may be indirect signatures of annihilating/decaying 
DM in our Universe. Although recent studies have claimed that the AMS-02 an¬ 

tiproton flux is within the uncertainties of the astrophysical secondary antiproton flux, 
the predicted secondary antiproton flux still tends to be smaller than the observed one at 
higher energy > 100 GeV, which may indicate a DM contribution in that energy range. 
Refs. [31 m |5] also derived upper bound on the DM signal, but the conservative bound |3] 
is still weak, which allows a large DM contribution. 

In this letter, we consider annihilating and decaying DM as a possible source of the 
AMS-02 antiproton flux. We show that annihilating/decaying DM with a mass of O(TeV) 
can explain the antiproton flux in the high energy range. We also investigate the impli¬ 
cations for supersymmetric (SUSY) DM. It is shown that the AMS-02 antiprotons may 
originate from Wino DM with a mass of 2-3 TeV. Surprisingly, the Wino mass of around 
3 TeV, which is suitable for the thermal relic DM scenario, can explain the observed an¬ 
tiproton data. Another interesting DM candidate is gravitino with an R-parity violation. 
It is shown that (9(TeV) gravitino DM with an R-parity violation can also be the source 
of the AMS-02 antiprotons. 

2 Antiproton from annihilating and decaying DM 

The flux o^rimary antiprotons from DM annihilation/decay at the Solar System, f = Uq, 
is given bjo 

( 1 ) 

where v{T) is the velocity of the antiproton with kinetic energy T and fp{T,f) is the 
antiproton number density per unit kinetic energy. The propagation of antiprotons is 
described by a cylindrical stationary diffusion model [6] 

= V [K{T)Vfp{T, f)] - ^ [sign{z)VJp{T, f)] - 2h5{z)T,^^{T)fp{T, f) + Q(T, f), (2) 

with boundary conditions fp(T, r) = 0 at r = i? and z = ±L, where (r, z) are the 
galactic cylindrical coordinates. Here, the effects of energy losses, reacceleration, and 
tertiary antiprotons are neglected. 

^We neglect the difference between the ffuxes at the top of the Earth’s atmosphere and at the 

interstellar <1>U due to the solar modulation, since its effect is 0(1)% for T > 50 GeV. 
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Model 

i?(kpc) 

T(kpc) 

iLo(kpcVMyr) 

(5 

K(km/s) 

MIN 

20 

1 

0.0016 

0.85 

13.5 

MED 

20 

4 

0.0112 

0.70 

12 

MAX 

20 

15 

0.0765 

0.46 

5 


Table 1: Propagation parameters [5]. 


In Eq. ([2]), K{T) is the diffusion coefficient and assumed to be spatially constant. 
It is parametrized as K{T,r) = K{T) = iPo/9(p/GeV)'^, where p and /5 = v{T)/c are 
the momentum and velocity of the antiproton, respectively. The Vc term represents the 
convective wind, which is assumed to be constant and perpendicular to the galactic plane. 
The third term represents the annihilation of the antiproton on interstellar protons in the 
galactic plane, where h represents the thickness of the galactic plane and rann(2") = 
(riH + 4^/^?7,He)c’‘p™'y(T) is the annihilation rate. We take h = 0.1 kpc, tt-h = 1 cm“^, 
UHe = O.OTuh, and a™ given in Refs. [3 [8], 


ann 

^PP 


661(1 + 0.0115T-°-^^^ - 0.948T° °i®i)mb T < IS.SGeV, 
36T-°-®mb T> IS.SGeV. 


(3) 


Lastly, Q{T,f) is the source term of the antiprotons. We adopt the set of propagation 
parameters R, L, Ko, 6, and W in Ref. [9], which are shown in Tabled! 

The source term for DM annihilation/decay is given by 


Q{T,f) = q(T)^^^ (4) 

where dNp{T)/dT is the energy spectrum of the antiproton per one annihilation/decay, 
and q{r) is given by 


q{r) 

q{r) 



PDMdfp 

moM 


1 f PDM(|f1) \ 

RDM \ / 


2 


for annihilating DM, 
for decaying DM. 


(5) 

( 6 ) 


Here, moM and poMd^l) are the mass and the density profile of the DM, respectively. In 
addition, (av) is the annihilation cross section for the annihilating DM case, while tdm is 
the DM lifetime for the decaying DM case. 

The differential equation ([2]) can be solved analytically, which leads to 



vjT) 

Air 


G{T) 


dNpiT) 

dT 


(7) 


For the source spectrum dNp{T) /dT, we consider the following DM annihilation and decay 
channels: 
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Figure 1: Energy spectra of antiproton produced by non-relativistic annihilation of DM. 
Red and blue lines correspond to those for yy; —)■ W~^W~ and bb, respectively. The DM 
mass is taken to be 2 TeV. 

• annihilation: yy —)> W~^W~, 

• decay: y —)■ 

where x denotes the DM and £ is a charged lepton. In Fig. [H we show the numerical 
result for the antiproton spectrum from yy —)■ W~^W~ obtained by Pythia 6.4 [TU]. Our 
results ^ree well with the htting formula in Ref. HU, which we use in the following 
analysiso The spectrum from the decaying DM with a mass of moM is given by that from 
the annihilating DM with a mass of mDM/2 rescaled by the factor of 1/2. 

In Fig. [U we also show the antiproton spectrum from yy —)■ bh. The source spectra for 
W~^W~ and bb are relatively close in the parameter range of our interest. We have checked 
that the resultant antiproton flux from yy —)■ bh is similar to the one from yy —)■ W~^W~. 
The analytic expression for G{T) is given by [6] 



( 8 ) 


where 






^ As for the fitting parameters Pi{muM) in Ref. [TT], we used PiirriDM = 5 TeV) for todm > 5 TeV. 
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Model 

Oo 

Cli 

0-2 

Os 

MIN 

1.1127 

1.7495 

-1.2730 

0.1412 

MED 

3.0662 

0.8814 

-0.8377 

0.09178 

MAX 

4.5815 

-0.3546 

-0.2322 

0.02524 


Table 2: fitting paramaters for decaying DM 


Here, Jq and Ji are the zeroth and first order Bessel functions of the hrst kind, respectively, 
and (i are the successive zeros of Jq. We have calculated the Green’s function G{T) by 
using the NFW density profile na 

\r\ \l + \r\/rsj 

with Pq = 0.3GeV/cni^, r© = 8.5 kpc and = 20kpc. We parametrize the result as 
~ 1 

G{T) = - (av) — ^G^^JT) annihilating DM, (10) 

2 

G{T) = — — -Gdec(T) decaying DM. (11) 

"^DM'TdM 

For annihilating DM, our numerical result of Gann(^) agrees well with the htting formula 
given in Ref. [ 13 ]. For decaying DM, our numerical result is well reproduced with the 
following htting function 

Gdec{T) = exp (oo + aiT + a 2 r^ + asr^) x lO^'^sec, (12) 

where 


r = logio(T/GeV), (13) 

and the coefhcients are shown in Table [2] for MIN, MED, and MAX propagation modelsH 
In Fig. [21 we show our numerical results of p/p ratio for annihilating and decaying 
DM. Here, for the proton hux, we adopt the following htting formula 


m-2sr-isec“^GeV“^ 


[10.0 — 9{—T3oq)3.0 Tsoo — ^(t3oo)0.6 Tsoo] x 10^ 



(14) 


where raoo = logio(2"/300 GeV), which reproduces the newly released proton hux by 
AMS-02 [ 2 ] well for T > 30GeV. In the hgures, we also show the background spectrum 
represented by the htting function in Ref. [T3] : 

f \ 

logic ^ ^ ,,-i = -1-64 + 0.07r - - 0.02 t ^ + 0.028rh (15) 

\ m“^sr“^sec“^GeV I 

^Our result slightly differs from the htting formula presented in Ref. m for high energy region, 
T > 0(100) GeV. 
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As can be seen from the figures, the antiproton flux can be explained by annihilat¬ 
ing/decaying DM with masses of (9(TeV). Notice that the cross sections and lifetimes 
used in Fig. [2] are not the best-fit values, but just for presentation. The signal antipro¬ 
ton flux is proportional to {av) or TjJjJj, and hence the AMS-02 antiproton flux can be 
explained in a wide range of DM mass if {av) or tdm is appropriately chosen. 


3 Implications for supersymmetric DM 

In this section, we briefly discuss the implications for some of SUSY DM candidates: 
annihilating Wino DM and decaying gravitino DM. 

3.1 Wino Dark Matter 

As shown in the previous section, if the antiproton flux is from the annihilating DM, it re¬ 
quires a mass of O(TeV) and relatively large annihilation cross section of 0(10“^^ cm^s“^). 
Such a parameter space is natural in the Wino DM scenario. In this case, the DM anni¬ 
hilation cross section is determined by its mass, and hence the antiproton flux depends 
only on the Wino mass. 

In Fig. [3], we show the p/p ratio for several Wino masses for MIN, MED and MAX 
propagation models. Here, we adopt the annihilation cross section in Ref. [T^ [TB] . As 
seen in the figure, Wino DM can account for the observed p/p ratio for several choices of 
masses. Interestingly, for the MED and MAX propagation models, Wino mass of 2.9 TeV 
can fit the antiproton data. Such a Wino mass is indeed predicted by the thermal relic 
abundance of Wino DM [min]. As can be seen in the figures, smaller Wino masses can 
also explain the data. In this case, a non-thermal production of Wino DM is necessary, 
such as gravitino [ISl [201 [211122], moduli [TSl 123], or Q-ball decay [23] • The antiproton 
flux at the AMS-02 may be the first hint of the Wino DM. 


3.2 Gravitino DM with R-parity violation 

In the scenario of gravitino DM with an R-parity violation [2S1 12S], a small R-parity 
violating coupling and the Planck-suppressed interaction lead to a long DM lifetime. In 
the case of bilinear R-parity violation, its lifetime is given by [221 126] 


r 3/2 ~ 10^® sec 


A \ ^ / ^3/2 
10-V VI TeV 



(16) 


where m 3/2 is the gravitino mass and A is the R-parity violating effective coupling. 

As shown in Fig. [2l if the high energy antiproton flux is from decaying gravitino DM, 
it implies 777 , 3/2 > 0{1) TeV and T 3/2 — 10^® — 10^^ s. This corresponds to an R-parity 
violating coupling of A ~ 10“^ —10“®. Such a value of R-parity violation is attractive from 
the cosmological point of view, since it is large enough to solve the constraint from the 
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Figure 2: The antiproton to proton ratio for MIN (top), MED (middle) and MAX (bot¬ 
tom) propagation models. The red lines are those predicted by the DM annihilation (left) 
and decay (right). For the annihilation case, the DM mass is taken to be 0.5 (solid), 
2 (dotted), 10 (dashed), and 20 TeV (long-dashed), while the annihilation cross section 
is taken to be 2 x 10“^^, 2 x 10“^^, and 6 x 10“^^ cm^/sec, for MIN, MED, and MAX 
propagation models, respectively. For the decay case, the DM mass is taken to be 1 
(solid), 3 (dotted), 10 (dashed), and 30 TeV (long-dashed), while the lifetime is 1 x 10^®, 
5 X 10^®, and 2 x 10^^ sec, for MIN, MED, and MAX propagation models, respectively. 
The background is shown in the green line, and the AMS-02 data are shown by the cyan 
points. 
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Figure 3: The antiproton to proton ratio in the Wino DM scenario for MIN (top), MED 
(middle) and MAX (bottom) propagation models. Red lines are those for the Wino mass 
of 2.9 TeV, while blue lines are those for the Wino mass of 2.2 TeV (MIN), 1.7 TeV (MED), 
and 1.2 TeV (MAX). The solid lines are signal plus background, while the dashed lines 
are signal-only. The background is shown in the green line, and the AMS-02 data are 
shown by the cyan points. 
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long-lived next-to-liglitest SUSY particle (NLSP) while small enough to avoid the baryon 
erasure in combination with sphaleron (cf. [27]). 

The gravitino abundance is given by [28l 123 


^ 3 / 2 / 1 ^ ~ 0 . 15(2 In 




/ "^ 3/2 \ ^ 

UTeVy VlO^GeVy ’ 


(17) 


where gs denotes the SU(3) gauge coupling constant, nig the gluino mass and Tr the 
reheating temperature after inflation. Therefore, for explaining the AMS-02 antiprotons 
with gravitino DM of m 3/2 > (a few) TeV, we need Tr <10® GeV. This is slightly lower 
than the temperature required by the standard thermal leptogenesis scenario [2111321IS3] , 
while non-thermal leptogenesis [31] can explain the observed baryon asymmetry. 


4 Discussion 

Now let us discuss various observational constraints on the annihilating/decaying DM 
scenario for explaining the AMS-02 antiprotons. 

• Cosmic microwave background (CMB): Recent Planck observation on the cosmic 
microwave background anisotropy constrains the DM annihilation cross section as 
/efr(cru)/mDM ^ 4 X lO”^® cm^s“^GeV“^ [35]) where fes ~ 0.3 (in the case of DM 
annihilation into W~^W~ or bb) is the effective fraction of the energy per DM anni¬ 
hilation that ionizes the hydrogen at the epoch of recombination [3S]- This leads to 
the following constraint 


W<lxlO-“cmV'(^). (18) 

For the MAX and MED model, this constraint is satished. For the MIN model, the 
required cross section is too large to satisfy this bound. This does not constrain the 
decaying DM model, because the energy injection around the recombination epoch 
is sufficiently small. Big-bang nucleosynthesis also constrains the annihilating DM 
model [371 fT 6 l 138] . but the constraint is weaker than that from GMB. 

• Gosmic-ray positron: The annihilating/decaying DM also yields high-energy cosmic- 
ray positrons. We have explicitly checked that the typical annihilating/decaying DM 
models explaining the AMS-02 antiprotons do not conflict with the PAMELA [33] 
and AMS-02 [iQ] positron measurements for the MED and MAX models. 

In the case of gravitino DM decaying into or a sizable positron flux is 

expected. In particular, when the gravitino decays into IF^e^, a sharp edge of the 
positron spectrum at the energy around 7713 / 2/2 may be an interesting signature of 
this scenario. 

^Note that the NLSP quickly decays with R-parity violation without producing gravitinos. 







• Gamma-rays (continuum): Fermi satellite searches for gamma-rays from dwarf spheroidal 
galaxies and puts severe constraint on the DM annihilation cross section and decay 
rate. The constraint reads {av) < (0.2 — 10) x 10“^^ cm^s“^ for DM mass of 1 — 10 TeV 

in the case of annihilating DM HU. For decaying DM, the diffuse gamma-ray back¬ 
ground gives stringent constraint [12] and the typical constraint reads tdm ^ 10^’^ s 
for DM mass of 1 — 10 TeV. For both annihilating and decaying DM cases to explain 
the AMS-02 antiproton data, the MED and MAX models can satisfy the constraint. 

• Gamma-rays (line): Annihilating or decaying DM also produces line gamma. For 
example, the Wino DM annihilates into Zy. For DM mass of C>(TeV), the HESS 
telescope gives stringent constraint on such line signals from the Galactic center in 
the case of DM annihilation [H]. However, the constraint signihcantly depends on 
the DM density profile. According to Ref. jH] , in which constraint from the gamma- 
ray line from Galactic center is derived for the case of Wino DM, MED and MAX 
parameters are allowed for mild coring of the DM density prohle around the Galactic 
center. 

In summary, we studied the cosmic-ray antiproton flux from DM annihilation and 
decay in light of the recent AMS-02 result. It is possible to explain the observed p/p 
ratio by adding DM contributions to the typical astrophysical background. In particular, 
we found that Wino DM with mass of around 3 TeV can successfully account for the 
antiproton data, which is consistent with the present DM abundance in the standard 
thermal freezeout scenario. Decaying gravitino DM heavier than a few TeV can also 
explain the data. 
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Note added 
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